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Hybrid bgb molecules in which the heavy bb pair is bound together by the excited gluon field 
g are studied using the Born-Oppenheimer expansion and quenched numerical simulations. The 
consistency of results from the two approaches reveals a simple and compelling physical picture for 
heavy hybrid states. 

PACS number(s): ll.15.Ha, 12.38. Gc, 12.39.Mk 



In addition to conventional hadrons, QCD predicts the 
existence of glueballs and hybrid states which contain ex- 
cited gluon fields. Hybrid mesons with heavy bb quark 
pairs are the most amenable to theoretical treatment. 
They are also experimentally accessible: early results 
from the CUSB and CLEO collaborations [f]j|] revealed 
a complex resonance structure between the BB threshold 
and 11.2 GeV in e+e" annihiliation, precisely where the 
lowest hybrid excitations are expected [|j. 

In this work, we determine the masses of the lowest 
bgb states. Heavy hybrid mesons can be studied not 
only directly by numerical simulation, but also using the 
Born-Oppenheimer expansion which is our primary guid- 
ance for the development of a simple physical picture. 
The Born-Oppenheimer picture was introduced for the 
description of heavy hybrid states in Refs. [Q,g and was 
applied using hybrid potentials first calculated in lattice 
QCD in Ref. Q. In this new study, we work to leading 
order in the expansion and neglect higher-order terms 
involving spin, relativistic, and retardation effects. We 
test the accuracy of the Born-Oppenheimer approach by 
comparison with high- precision results from simulations. 

Our hybrid meson simulations are the first to exploit 
anisotropic lattices with improved actions; preliminary 
reports on some of our results have appeared previously 
W . The hybrid meson mass uncertainties with improved 
anisotropic lattice technology are dramatically smaller 
than those obtained in recent isotropic lattice studies in 
the nonrelativistic formulation of lattice QCD (NRQCD) 
using the Wilson gauge action Wm. We report here our 
final analysis on four distinct hybrid bgb states. Al- 
though the effects of dynamical sea quarks are not in- 
cluded in our quenched simulations, we will comment 
on their impact on the hybrid spectrum. The mass of 
the lowest hybrid cgc state was determined recently Q 
without NRQCD expansion for the slowly moving heavy 
quark and agrees with our Born-Oppenheimer results 0] 
(see caption of Fig. |l|) . 

The hybrid meson can be treated analogous to a di- 
atomic molecule: the slow heavy quarks correspond to 



the nuclei and the fast gluon field corresponds to the 
electrons Q. First, one treats the quark Q and anti- 
quark Q as spatially-fixed color sources and determines 
the energy levels of the excited gluon field as a function of 
the QQ separation r; each of these excited energy levels 
defines an adiabatic potential Vq q(i')- The quark mo- 
tion is then restored by solving the Schrodinger equation 
in each of these potentials. Conventional quarkonia are 
based on the lowest-lying static potential; hybrid quarko- 
nium states emerge from the excited potentials. Once 
the static potentials have been determined (via lattice 
simulations) , it is a simple matter to determine the com- 
plete spectrum of conventional and hybrid quarkonium 
states in the leading Born-Oppenheimer (LBO) approxi- 
mation. This is a distinct advantage over meson simula- 
tions which yield only the very lowest-lying states, often 
with large statistical uncertainties. In addition, the LBO 
wave functions yield valuable information concerning the 
structures and sizes of these states which should greatly 
facilitate phenomenological applications. 

The energy spectrum of the excited gluon field in the 
presence of a static quark-antiquark pair has been deter- 
mined in previous lattice studies 0. The three lowest- 
lying levels are shown in Fig. n^. These levels correspond 
to energy eigenstates of the excited gluon field charac- 
terized by the magnitude A of the projection of the to- 
tal angular momentum Jg of the gluon field onto the 
molecular axis, and by 77 = ±1, the symmetry quan- 
tum number under the combined operations of charge 
conjugation and spatial inversion about the midpoint 
between the quark and antiquark of the QgQ system. 
Following notation from molecular spectroscopy, states 
with A = 0, 1, 2, . . . arc typically denoted by the capi- 
tal Greek letters S,H, A, . . ., respectively. States which 
are even (odd) under the above-mentioned parity-charge- 
conjugation operation are denoted by the subscripts g 
(u). There is an additional label for the S states; S states 
which are even (odd) under a reflection in a plane con- 
taining the molecular axis are denoted by a superscript -I- 
(— ). In Ref. 0], the potentials are calculated in terms of 




FIG. 1. (a) Static potentials and radial probability densi- 
ties against quark-antiquark separation r for r;^^ — 450 MeV. 
The Eg potential becomes the familiar running Coulomb law 
as r becomes very small, (b) Spin-averaged bb spectrum in the 
LBO approximation (light quarks neglected) . Solid lines indi- 
cate experimental measurements. Short dashed lines indicate 
the S and P state masses obtained from the Schrodinger equa- 
tion with the S^ potential for Mb — 4.58 GeV. Dashed-dotted 
lines indicate the hybrid quarkonium states obtained from the 
IIu (L = l,2, 3) and "E^ (L = 0, 1, 2) potentials. Repeating 
the same analysis in the cgc system, we find the lowest IIu 
hybrid 1.19 GeV above the spin averaged ground state. 

the hadronic scale parameter tq |11|; the curves in Fig. n^ 
assume r^^ = 450 MeV (see below). Note that as r be- 
comes small (below 0.1 fm), the gaps between the excited 
levels and the S+ ground state will eventually exceed the 
mass of the lightest glueball. When this happens, the ex- 
cited levels will become unstable against glueball decay. 

Given these static potentials, the LBO spectrum is eas- 
ily obtained by solving the radial Schrodinger equation 
with a centrifugal factor (L'L ) = L(L + 1) - 2A2 + (Jl), 
where Lqq is the orbital angular momentum of the 
quark-antiquark pair. For the E+ potential, (Jg) = 0. 
For the IIu and E~ levels, we attribute the lowest non- 
vanishing value (Jg) = 2 to the excited gluon field. Let 
S be the sum of the spins of the quark and antiquark, 
then the total angular momentum of a meson is given by 
J = L + S. In the LBO approximation, the eigenvalues 
L(L + 1) and S(S -I- 1) of L^ and S^ are good quantum 
numbers. The parity P and charge conjugation C of each 
meson is given in terms of L and S by P = e ^_i)L+a+i 
and C = e ?7 (-1)^+^+8^ ^here L > A and e = 1 for 11+, 
e = — 1 for S^, and e = ±1 for A > 0. Note that for each 
static potential, the LBO energies depend only on L and 
the radial quantum number n. 

Results for the LBO spectrum of conventional bb and 
hybrid bgb states are shown in Fig. y. The heavy quark 
mass Mb is tuned to reproduce the experimentally-known 
T(1S) mass: Mx = 2Mb + Eo, where Eq is the energy of 
the lowest-lying state in the E+ potential. Level split- 
tings are insensitive to small changes in the heavy quark 



mass. For example, a 5% change in Mb results in changes 
to the splittings (with respect to the IS* state) ranging 
from 0.1 - 0.8%. 

Below the BB threshold, the LBO results are in very 
good agreement with the spin-averaged experimental 
measurements of bottomonium states. Above the thresh- 
old, agreement with experiment is lost, suggesting signif- 
icant corrections from higher order effects and possible 
mixings between the states from different adiabatic po- 
tentials. The mass of the lowest-lying hybrid (from the 
IIu potential) is about 10.9 GeV. Hybrid mesons from 
all other hybrid potentials are significantly higher lying. 
The radial probability densities for the conventional IS 
and IP states are compared with that of the lowest-lying 
IIu hybrid state in Fig. ^ Note that the size of the hy- 
brid state is large in comparison with the IS and IP 
states. For all of the hybrid states studied here, the 
wave functions are strongly suppressed near the origin 
so that the hybrid masses cannot be affected noticeably 
by the small-r instability of the excited-state potentials 
from bgb ^ bb -|- glueball decay. 

The applicability of the Born-Oppenheimer approx- 
imation relies on the smallncss of retardation effects. 
The difference between the leading Born-Oppenhcimer 
Hamiltonian and the lowest order NRQCD Hamiltonian 
is the p • A coupling between the quark color charge in 
motion and the gluon field. This retardation effect, which 
is not included in the LBO spectrum, can be tested by 
comparing the LBO mass splittings with those deter- 
mined from meson simulations in NRQCD. 

In order to obtain the masses of the first few excited 
hybrid states in a given symmetry channel, we obtained 
Monte Carlo estimates for a matrix of hybrid meson cor- 
relation functions Cij(t) = (0|Mi(t)M|(0)|0) at two dif- 
ferent lattice spacings. Because the masses of the hybrid 
mesons are expected to be rather high and the statistical 
fluctuations large, it is crucial to use anisotropic lattices 
in which the temporal lattice spacing at is much smaller 
than the spatial lattice spacing ag. Such lattices have 
already been used to dramatically improve our knowl- 
edge of the Yang-Mills glueball spectrum ||l2[. In our 
simulations, the gluons are described by the improved 
gauge-field action of Ref. [Q. The couplings /3, input 
aspect ratios (,, and lattice sizes for each simulation are 
listed in Table ||. Following Ref. ||l2[ , we set the mean 
temporal link ut = 1 and obtain the mean spatial link Ug 
from the spatial plaquette. The values for ro in terms 
of as corresponding to each simulation were determined 
in separate simulations. Further details concerning the 
calculation of ro/ag are given in Ref. p3]. Note that we 
set the aspect ratio using ag/at = ^ in all of our calcula- 
tions. By extracting the static-quark potential from Wil- 
son loops in various orientations on the lattice |l3], we 
have verified that radiative corrections to the anisotropy 
ag/at are small. The heavy quarks are treated within 
the NRQCD framework iQ, modified for an anisotropic 
lattice. The NRQCD action includes only a covariant 



TABLE I. Simulation parameters and results. The second 
errors listed in the results in the bottom six rows are due to 
uncertainties in setting the heavy quark mass. 



(/3,C) 




(3.0,3) 


(2.6,3) 


u* 




0.500 


0.451 


lattice 




15^x45 


10^ X 30 


# config! 


5, sources 


201, 16080 


355, 17040 


ro/as 




4.130(24) 


2.493(9) 


(C,nc) 




(0.25, 15) 


(0.15, 10) 


as Mo 




2.56 


3.90 


asM^i^ 




5.03(2) 


8.21(1) 


ro 5(1P - 


-IS) 


0.959(8)(3) 


0.998(6)(3) 


ro 5(2S - 


'IS) 


1.303(11)(10) 


1.252(8)(10) 


ro 5(Hi - 


-IS) 


3.287(53)(20) 


3.338(54)(20) 


ro (5(H2 - 


-IS) 


3.37(13)(1) 


3.443(47)(10) 


ro 5(H3 - 


-IS) 


4.018(55)(12) 


4.034(76)(12) 


ro S{R[ ^ 


-IS) 


4.204(67)(21) 


4.229(62)(21) 



temporal derivative and the leading kinetic energy op- 
erator (with two other operators to remove 0(at) and 
0(ag) errors); relativistic corrections depending on spin, 
the chromoelectric E and chromomagnetic B fields, and 
higher derivatives are not included. 

Our meson operators Mi(t) are constructed on a given 
time-slice as follows. First, the spatial link variables are 
smeared using the algorithm of Ref. 11^] in which ev- 
ery spatial link Uj(x) on the lattice is replaced by it- 
self plus ( times the sum of its four neighboring spatial 
staples, projected back into SU(3); this procedure is it- 
erated n^; times, and we denote the final smeared link 
variables by Uj(x). Next, let tpix) and xW denote the 
Pauli spinor fields which annihilate a heavy quark and 
antiquark, respectively. Note that the antiquark field is 
defined such that C'(/'(x)C^ = icryX*(x), where C is the 
charge conjugation operator. We define a smeared quark 

field by V'(x) = (l + g a^ A*^^) J ■0(x), (and similarly for 

the antiquark field) where g and n^ are tunable param- 
eters (we used g = 0.12,0.14 and n^ = 2 — 7) and the 
covariant derivative operators are defined in terms of the 
smeared link variables Uj(x). These field operators, in 
addition to the chromomagnetic field, are then used to 
construct our meson operators, which are listed in Ta- 
ble O. The standard clover-leaf definition of the chro- 
momagnetic field B is used, defined also in terms of the 
smeared link variables. Note that four operators are used 
in each of the ^ and 1 sectors. Because our NRQCD 
action includes no spin interactions, we use only spin- 
singlet operators. We can easily couple these operators 
to the quark-antiquark spin to obtain various spin-triplet 
operators, and the masses of such states will be degener- 
ate with those from the spin-singlet operators. 

In each simulation, the bare quark mass agMo is set by 
matching the ratio R = M^i„/(S(1P - IS), where M^j^^ is 
the so-called kinetic mass of the IS state and S{1P — IS) 



TABLE IL The meson spin-singlet operators used in each 
total angular momentum J, parity P, and charge conjugation 
C channel. Note that p = 0, 1, 2, and 3 were used to produce 
four distinct operators in the '' and 1~~ sectors. In the 
third column are listed the spin-triplet states which can be 
formed from the operators in the last column; the states in 
each row are degenerate for the NRQCD action used here. 



J"" 




Degeneracies 


Operator 


0-+ 


S wave 


1"- 


Xt At^)- " 4, 


1+- 


P wave 


0++,l++,2++ 


y} Li, 


1-- 


Hi hybrid 


0-+,l-+,2-+ 


xtB[A(2)]'' vi 


1++ 


H2 hybrid 


0+-,l+-,2+- 


X^ BxA i/) 


0++ 


H3 hybrid 


1+- 


X^ B-A ■ip 



is the energy separation between the IS and IP states, 
to its observed value 21.01(6). The kinetic mass M^j^^ is 
determined by measuring the energy of the IS state for 
momenta p= (0,0,0), 27r(l,0, 0)/L, and 27r(l, 1, 0)/L, 
where L is the spatial extent of the periodic lattice. These 
three energies are then fit using Eq -I- p^/(2M^;jj) to ex- 
tract M^jjj. Several low statistics runs using a range of 
quark masses were done in order to tune the quark mass. 
From the results of these runs, we estimate that the un- 
certainty in tuning the quark mass is about 5%. 

The simulation results are listed in Table ffl The 
masses mj in the 1^ , l"*"^, and 0+"*" channels are ex- 
tracted by fitting the single correlators Ci(t) to their ex- 
pected asymptotic form Ci(t) -^ Ziexp(— niit) for suffi- 
ciently large t. In each of the '" and 1 channels, 
we obtain a 4 x 4 correlation matrix. The variational 
method is then applied to reduce the 4x4 matrix down 
to an optimized 2x2 correlation matrix C°j^ (t). For suf- 
ficiently large t, we fit all elements of this matrix using 
^ QZipZjpexp(— mpt) to extract the two lowest-lying 
masses. In this way, we obtain an estimate of the 2S 
mass, as well as the first excited hybrid state H'j^. The 
effective masses corresponding to several of the correla- 
tion functions obtained in the (3 = 3.0, ^ = 3 simulation 
are shown in Fig. 0. 

The simulation results for the level splittings (in terms 
of ro and with respect to the IS state) are shown in Fig. 
against the lattice spacing. Small finite-as errors are ev- 
ident in the IP and 2S splittings from the coarse lat- 
tice simulation; none of the four hybrid splittings show 
any significant discretization errors. The simulation re- 
sults compare remarkably well with the LBO predictions, 
shown as horizontal lines in Fig. 0. In the LBO approxi- 
mation, the Hi and H2 mesons correspond to degenerate 
lPn„ states of opposite e, the H3 hybrid corresponds to 
a ISj,- state, and the H'^ corresponds to a 2Pnu level; 
furthermore, the II3 and H'^ hybrids are predicted to be 
nearly degenerate, with the H'^ lying slightly lower. The 
simulation results share these same qualitative features, 
except that the H'^ lies slightly higher than the II3. The 
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FIG. 2. Effective masses iiicff(t) 
for two of the hybrids from the (/3, ^) - 



= -ln[C(t + at)/C(t)] 
(3.0, 3) simulation. 



LBO approximation reproduces all of the level splittings 
to within 10%. In Fig. g, we also show results ||l^ for 
the IP and 2S splittings for an NRQCD action includ- 
ing higher order relativistic and spin interactions; the 
effects of such terms are seen to be very small. Note that 
spin-dependent mass splittings are difficult to estimate in 
hybrid states since the excited gluon field extends on the 
scale of the confinement radius with a nonperturbative 
wavefunction when its color magnetic moment interacts 
with the heavy-quark spins. 

To convert our mass splittings into physical units, we 
must specify the value of rp. Using the observed value 
for the IP - IS splitting, we find that r^"^ = 467(4) MeV; 
using the 2S — IS splitting, we obtain r^^ = 435(5) MeV. 
This discrepancy is caused by our neglect of light quark 
effects 0- Taking r^^ == 450(15) MeV, our lowest-lying 
hybrid state lies 1.49(2)(5) GeV (the second error is the 
uncertainty from ro) above the IS state. 

Hybrid and conventional states substantially extending 
over 1 fm in diameter are vulnerable to light-quark vac- 
uum polarization loops which will dramatically change 
the static potentials through configuration mixing with 
BB mesons; instead of rising indefinitely with r, these po- 
tentials will eventually level off since the heavy QgQ state 
can undergo fission into two separate Qq color singlets, 
where q is a light quark. We expect that the plethora of 
hybrid and conventional states above 11 GeV obtained 
from the quenched potentials will not survive this split- 
ting mechanism as observable resonances. For quark- 
antiquark separations below 1.2 fm or so, there is evi- 
dence from recent studies p8|,p| that light-quark vacuum 
polarization effects do not appreciably alter the E^ and 
IIu inter-quark potentials (for light-quark masses such 
that uItt ^ nip/2)- Since such distances are the most 
relevant for forming the lowest-lying bound states, the 
survival of the lightest bgb hybrids as well-defined reso- 
nances above the BB threshold remains conceivable. 

During the preparation of this work we learned about 
new results [ |l9[ which have considerable overlap with our 
NRQCD simulations. 
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FIG. 3. Simulation results for the level splittings (in terms 
of ro and with respect to the IS state) against the lattice 
spacing as. Results from Ref. [[I6| using an NRQCD action 
with higher-order corrections are shown as □ and A; all other 
symbols indicate results from this work. Some points have 
been shifted horizontally to prevent overlaps. The horizontal 
lines show the LBO predictions. 
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